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Abstract 

The energetics of the gas-phase SAr2 reactions Y~ + CH3X — > CH3Y + 
X~ (where X,Y = F, CI, Br), were studied using (variants on) the recent 
Wl and W2 ab initio computational thermochemistry methods. These cab 
culations involve CCSD and CCSD(T) coupled cluster methods, basis sets of 
up to spdfgh quality, extrapolations to the one-particle basis set limit, and 
contributions of inner-shell correlation, scalar relativistic effects, and (where 
relevant) first-order spin-orbit coupling. Our computational predictions are 
in excellent agreement with experimental data where these have small error 
bars; in a number of other instances re-examination of the experimental data 
may be in order. Our computed benchmark data (including cases for which 
experimental data are unavailable altogether) are used to assess the quality of 
a number of semiempirical compound thermochemistry schemes such as G2 
theory, G3 theory, and CBS-QB3, as well as a variety of density functional 
theory methods. Upon applying some modifications to the level of theory 
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used for the reference geometry (adding diffuse functions, replacing B3LYP 
by the very recently proposed mPWlK functional [Lynch, B.J.; Fast, P.L.; 
Harris, M.; Truhlar, D.G. J. Phys. Chem. A 2000, IO4, 4811]), the com- 
pound methods appear to perform well. Only the 'half-and-half functionals 
BH&HLYP and mPWH&HPWQl, and the empirical mPWlK functional, con- 
sistently find all required stationary points; the other functionals fail to find 
a transition state in the F/Br case. BH&HLYP and mPWH&HPW91 some- 
what overcorrect for the tendency of B3LYP (and, to a somewhat lesser extent, 
mPWlPWQl) to underestimate barrier heights. The Becke97 and Becke97-1 
functionals perform similarly to B3LYP for the problem under study, while the 
HCTH and HCTH-120 functionals both appear to underestimate central bar- 
riers. HCTH underestimates complexation energies; this problem is resolved 
in HCTH-120. mPWlK appears to exhibit the best performance of the func- 
tionals considered, although its energetics are still inferior to the compound 
thermochemistry methods. mPWlK, however, appears to be very suitable for 
generating reference geometries for more elaborate thermochemical methods 
in kinetics applications. 
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I. INTRODUCTION 



Due to the central importance of bimolecular nucleophilic substitution (SAf2) reactions 
in organic chemistry,0i the prototype S Ar2 reactions 

Y- + CH3X — > CH3Y + X" (X, Y = F, CI, and Br) (1) 

have aroused considerable interest in the past three decades. (Halomethanes have also 
received considerable attention in the area of atmospheric chemistry in connection with 
global warmingil and ozone layer destructionE) Both theoretical and experimental studies 
see Refs.ii for reviews) indicate that the preferred gas phase reaction pathway involves a 
backside attack of the halide ion, Y^, at the carbon atom followed by the familiar 'Walden 
inversion' of the CH3 group. The resulting reaction profile (Figure 1) exhibits two local 
minima, i.e. entry and exit channel ion-molecule complexes Y"- ■ -CHsX and YCHs- ■ -X", 
connected by a central transition state [Y- ■ -CHa- ■ -X]^, which has D^h and C^^ symmetries 
in the identity (X=Y) and nonidentity (Xt^Y) cases, respectively. Although the qualitative 
form of this reaction pathway is widely accepted for substitution reactions in the gas phase, 
there is still considerable uncertainty about the exact energetics. 

Bohme et a/.i and Brauman et a/.i were the first to investigate the gas-phase SAr2 re- 
actions experimentally. Brauman and coworkers concluded that the measurements were 
best explained by a double-well potential with a central barrier. Subsequent experimental 
studiesiHli for a series of anionic nucleophiles with alkyl halides revealed that changes in 
the nucleophile, leaving group, and alkyl moiety leads to a wide variation of reaction rate 
constants; their observed variation was attributed to the central barrier height. 

The double-well Sjv2 potential energy surface also finds abundant theoretical support 
from ab initio calculations, which are currently one of the most useful tools for evaluating 
reaction potential energy profiles. For instance, Chandrasekhar et al^ presented a compre- 
hensive examination of Cl/Cl identity SAr2 reaction at the HF/6-31G* level, and Tucker and 
TruhlaiO examined the SAr2 reactions at the MP2/6-31G* level. Wladkowski et al.EJ studied 
the F/F identity Sjv2 reaction by large-scale coupled cluster theory involving single and dou- 
ble excitation operators with an a posteriori a quasiperturbative treatment of the effects of 



connected triple substitutions (CCSD(T)) (see Refs.ll3llZl for reviews). Later, Radom, Pross 
and coworkers have carried out ab initio molecular orbital calculations at the G2(+) level 
of theory for the back-side identityS and nonidentityEl SAr2 reactions. These authors have 
also investigated the identity front-side SAr2 reactions with retention of configuration.^ The 
G2(+) theory is essentially G2 theory carried out from MP2/6-31+G* (rather than MP2/6- 
31G*) geometries and employing scaled HF/6-31+G* (rather than HF/6-31G*) zero-point 
energies. For the bromine and iodine containing systems, these authors employed Hay- 
WadtS relativistic effective core potentials (RECPs). Botschwina and coworkers examined 
the stationary points of the potential surface for the F/Cl nonidentity SAr2 reactionil and 
for the Cl/Gl identity Sj\f2 reactionll by means of large-scale CGSD(T) calculations. Finally, 
a referee brought a very recent large-scale coupled cluster study by Schmatz et al@ to our 
attention. 

Despite the well known successes (e.g. of the increasingly popular DFT (density 
functional theory) methods,^ their performance for transition state structures and reac- 
tion barrier heights leaves something to be desired. For instance, DurantJ^l found that the 
B3LYP, B3P86 and B3PW91 functionals all systematically underestimated barrier heights, 
while only the Becke half-and-half/Lee, Yang and Parr (BH&HLYP)ii functional predicted 
transition state barrier heights reasonably well — despite the fact that its performance for 
thermochemical and other properties is generally substantially poorer than that of B3LYP 
and B3PW91. Baker et al^ arrived at a similar conclusion stating that the currently avail- 
able density functionals are unable to provide a correct description of the transition states. 
For the prototype 8^2 reactions (Gl/Cl and Gl/Br) considered here, Radom and coworkers^ 
found that the popular B3LYP@'i exchange-correlation functional significantly underes- 
timated the overall and central barrier heights compared to the G2(+) and experimental 
results. 

Nevertheless, the size of the systems involved in kinetic and mechanistic problems of 
organic and organometallic interest often makes DFT the only practical option. As a mat- 
ter of fact, our group has recently reported DFT studies of the mechanism of competitive 
intramolecular C-G and G-H bond activation in rhodium(I) pincer complexes^ and of the 
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Heck reaction.E3 

Aside from BH&HLYP, better performance for barrier heights has been claimed for a 
number of newer exchange-correlation functionals. For example, Adamo and BaroneS found 
that their mPWlPW91 (modified Perdew-Wang 1991 1-parameter hybrid exchange with 
Perdew-Wang 1991 correlation^) at least correctly predicts a positive overall barrier for the 
Cl/Cl identity SAr2 reaction, although it is still being underestimated. Very recently, Truhlar 
and coworkersii proposed a new hybrid model called the modified Perdew-Wang 1-parameter 
model for kinetics (mPWlK). In this empirical functional, the coefficient X for admixture 
of "exact" Hartree-Fock exchange 



(where X=l/4 for standard mPWlPW91) was determined (using the fairly small 6-31+G* 
basis set) by minimizing the average deviation from a set of 40 barrier heights (20 forward, 
20 reverse) obtained from a combination of experiment and theory (see Ref.S for details). 
(Note that the Walden inversion, or for that matter cationic or anionic reactions of any 
kind, were not part of the parametrization set.) It was found that mPWlK reduced the 
mean unsigned error in reaction barrier heights by a factor of 2.4 over mPWlPW91 and by 
a factor of 3 over B3LYP. 



Kassel-Marcus (RRKM)EI theory were also employed to examine the Siv2 reactions. Results 
from such studies (se and references therein) suggested that the assumption of statistical 
behavior in ion-molecule intermediate complexes is not valid. This "nonstatistical" behavior 
has been documented for several halide-methyl halide reactions and a thorough discussion 
is given by Hase.il Classical trajectory simulations performed by Hase and coworkers^ ques- 
tioned the basic assumptions of statistical theories and found that the trajectory calculations 
are very useful in interpreting the kinetics and dynamics of S Ar2 reactions. 

Despite the enormous amount of work in the past, there are still significant gaps in 
the experimental data for the gas-phase SAr2 reactions and, even where data are available, 
the results often possess large uncertainties. Recently, two computational thermochemistry 
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methods known as Wl and W2 (Weizmann-1 and Weizmann-2) theorycj have been developed 
in our laboratory. These are free of parameters derived from experiment and on average 
can claim 'benchmark accuracy' (defined in Ref.0 as a mean absolute error of 1 kJ/mol, 
or 0.25 kcal/mol) for molecular total atomization energies (TAEs) of first-and second-row 
compounds. The primary objective of the present study is to obtain high-quality energetic 
data for reaction (1) by means of Wl and W2 theory. Using these benchmark data, we 
shall then examine the performance of various DFT methods and ab initio computational 
thermochemistry methods such as Gl,@ G2,§ G3,0 and CBS-QB3§ theories. 



II. COMPUTATIONAL METHODS 

All calculations were carried out on the 4-processor Compaq ES40 of our research group, 
and on the 12-CPU SGI Origin 2000 of the Faculty of Chemistry. 

Energetics for the gas-phase stationary points for all six surfaces (i.e. F/F, Cl/Cl, Br/Br, 
F/Cl, F/Br, and CI/Br) were obtained by means of the Wl' method described in Refs.SS. 
The Wl' method@ is a minor variation on Wl theory^ that exhibits improved accuracy 
for second-row systems at no additional computational cost. For a detailed description and 
theoretical and empirical arguments for each step, see RefEi; we shall merely summarize 
the main points here for the sake of clarity. The basis sets employed are mostly Dunning's 
augmented correlation consistent ra-tuple zeta&^ (aug-cc-pVnZ) basis sets; for second-row 
atoms high-exponent d and / functions were added (denoted '+2d' or '+2dlf') as recom- 
mended in Ref.ii for accommodating inner polarization. Since the standard aug-cc-pVnZ 
basis sets for bromine!! already contain quite high-exponent d functions in order to describe 
the 3d orbitals, no 'inner polarization' functions were deemed to be necessary on Br. We 
may distinguish the following six components in the 'bottom-of-the-well' TAE at the Wl' 
level: 

• The SCF component of the TAE is obtained using the aug-cc-pVDZ+2d, aug-cc- 
pVTZ+2d, and aug-cc-pVQZ+2dlf basis sets, and extrapolated to the infinite-basis 
limit using the geometric expressionii A + B ■ C~^, where the 'cardinal number' 
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L={2,3,4} for these three basis sets. (It is identical to the maximum angular mo- 
mentum present for nonhydrogen atoms. Regular cc-pVnZ basis sets were used on 
hydrogen atoms throughout.) 

• The CCSD (coupled cluster with all singles and double^!) valence correlation contri- 
bution to TAE is obtained using the aug-cc-pVTZ+2d and aug-cc-pVQZ+2dlf basis 
sets, then extrapolated to the infinite basis hmit using the expression A + B / L^ "^"^ . 

• The (T) connected triple excitations componentil of TAE was computed using the aug- 
cc-pVDZ+2d and aug-cc-pVTZ+2d basis sets, and extrapolated to the infinite basis 
limit using the expression A + B /L^'^^. 

• The inner-shell correlation contribution was computed as the difference between 
CCSD(T)/MTsmall0 values with and without constraining the inner-shell orbitals 
to be doubly occupied. The very deep-lying chlorine (Is) and bromine (ls,2s,2p) or- 
bitals were doubly occupied throughout; the 'inner-shell correlation' thus represents 
carbon (Is), chlorine (2s, 2p) and bromine (3s,3p,3d) correlation. (Basis set superposi- 
tion error, BSSE, can be an issue for inner- shell correlation energies in heavier element 
systems;0 our experienced! suggests that BSSE in the W1/W2 inner shell correlation 
contributions largely cancels with basis set incompleteness.) 

• The scalar relativistic contribution was computed as expectation values of the one- 
electron Darwin and mass- velocity (DMV) operator J^il for the ACPF/MTsmall (av- 
eraged coupled pair functionalii) wave function, with all inner-shell electrons correlated 
except for chlorine (Is) and bromine (ls,2s,2p). 

• The spin-orbit contribution to TAE, in the present case of all-closed-shell systems, is 
nothing more than the sum of the atomic fine structure corrections. 

Where our computational hardware permitted (in practice, for F/F, Cl/Cl, and the 
Br/Br transition state), we also carried out even more demanding W2h calculations. In W2 
theory, the same steps occur as above, except that the three valence basis sets are aug-cc- 
pVTZ+2dlf, aug-cc-pVQZ+2dlf, and aug-cc-pV5Z+2dlf (with L=3,4, and 5, respectively) 
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and that the extrapolation formulaic used for the CCSD and (T) steps is simply A + B/L^. 
The W2h variantH indicates, in this particular case, the use of unaugmented cc-pVnZ basis 
sets on carbon. The largest basis set CCSD step was carried out using the direct algorithm of 
Lindh, Schiitz, and Werner.S All these calculations were performed using MOLPRO 98. 1@ 
and a driver for the W1/W2 calculation^! written in MOLPRO's scripting language. Ref- 
erence geometries were obtained primarily using the BSLYpHii density functional method, 
which employs the Lee-Yang-Parr0 correlation functional in conjugation with a hybrid ex- 



change functional first proposed by Becke.E^ 

A number of lower-level procedures were validated against the Wl' and W2h re- 
sults. These include the following set of DFT exchange- correlation functionals and ba- 
sis sets: B3LYP/cc-pVTZ(+X), BH&HLYP/cc-pVTZ(+X), mPWlPW91/cc-pVTZ(+X), 
mPWlK/6-31+G*, mPWlK/cc-pVDZ(+X), mPWlK/cc-pVTZ(+X), where (+X) indicates 
that diffuse functions are included only for halogens. BH&HLYPi is essentially the B3LYP 
method, with the exception that the fraction of HF exchange is 50% (H&H denotes "half and 
half). Analogous to BH&HLYP, we have also performed mPWH&HPW91.i^ In addition, we 
carried out calculations using the standard Gl, G2, G3, and CBS-QB3 model chemistries. 
The MP2/6-31G* and B3LYP/6-311G(2d,d,p) levels of theory (used for the reference ge- 
ometries in G3@ and CBS-QB3,0 respectively) fail to find stationary points for several of 
the ion-molecule complexes in the nonidentity cases (because of the absence of diffuse func- 
tions). Therefore, we have defined, by analogy with Radom and coworkersjll^ G3(-|-) and 
CBS-QB3(+) model chemistries where MP2/6-31+G* and B3LYP/6-311+G(2d,d,p) refer- 
ence geometries, respectively, are used. The G2(+) results quoted in the tables are taken 
from Radom and coworker s.llilli All of these calculations were carried out using Gaussian 98 
rev. or trivial modifications thereof. Following the recommendations in Ref.ii, larger 

grids than the default were used in the DFT calculations if necessary, specifically a pruned 
(99,590) grid for integration and gradients, and a pruned (50,194) grid for the solution of 
the coupled perturbed Kohn-Sham equations. 

Finally, following very recent suggestions in the literature^ that some of these function- 
als may perform better for transition states, some calculations using the novel B97 (Becke- 
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1997) JZ2l B97-1 (reparametrized Becke-1997),0 HCTH (Hamprecht-Cohen-Tozer-Handy),0 
and HCTH-120 (reparametrization of HCTH including anions and weakly interacting 
systems)^ functionals were carried out by means of a slightly modified version of NWCHEM 
3.3.1.il 

III. RESULTS AND DISCUSSION 

A. Reference Geometries 

Reference geometries for the Wl' calculations were mostly obtained at the B3LYP/cc- 
pVTZ+1 level, where the '+1' signifies the addition of a high-exponent d function on second- 
row elements.0 (The Br basis set already includes high-exponent d functions to cover the (Be?) 
orbital.) It was previously shownil that B3LYP/cc-pVTZ geometries for stable molecules 
are generally within a few thousandths of an A from experiment, as well that the use of 
B3LYP/cc-pVTZ+l rather than much costlier CCSD(T)/cc-pVQZ+l reference geometries 
insignificantly affects computed energies.El (Modifications of popular computational ther- 
mochemistry methods that use DFT reference geometries include variant of G2 theory, 
G3//B3LYP,0 and CBS-QB3.0) In some cases where B3LYP fails to locate the required 
stationary point, we used mPWlK/cc-pVTZ(+X) reference geometries. For the W2h cal- 
culations, CCSD(T)/cc-pVQZ+l reference geometries were used. The geometries of all the 
structures involved in the present study calculated at various levels of theory are provided 
in the Supplementary Material. 

B. Energetics 

In order to assess the accuracy of Wl' and W2h results, we consider first the total 
atomization energies (TAEs) of CH3X and electron affinities (EAs) of X~. A summary of 
our computed results and their components for the react ants/pro ducts of the S7v2 reactions is 
presented in Table 1. The final energies presented in the last column of the Table correspond 
to EAs of X (X = F, CI and Br) and TAEs without zero-point vibrational energy (ZPVE) 
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of CH3X. The inner-shell correlation contributions are all positive and the largest is 1.49 
kcal/mol for CHsBr. The core correlation contribution for TAE(CH3X) is found to increase 
in the order F < CI < Br, while for EA(X) it increases in the order CI < F < Br. The 
importance of Darwin and mass-velocity corrections increases, as expected, with increasing 
atomic number (Z) of X and its contribution becomes substantial when X = Br. 

It is perhaps more pertinent for our purposes to examine the relative energies (with 
respect to reactants) of the ion-molecule complexes and transition state structures. Table 
2 presents Wl' and W2h results for identity reactions and only Wl' results for nonidentity 
reactions. W2h calculations for the nonidentity reactions are extremely expensive as the 
reaction intermediates are less symmetric. Moreover, the size of the bromine atoms prevents 
us from performing a W2h calculation on the identity Br^- ■ -CHsBr ion-molecule complex. 
Likewise, we could not obtain the core correlation contributions for the Br^- ■ -CHsBr ion- 
molecule complex at the Wl' level of theory. It was previously establishecS that the inclusion 
of connected triple excitations in CCSD(T) is absolutely necessary for reliable core correlation 
contributions: the n^N'^ CPU time dependence of the (T) step dominates the required CPU 
time for CI and Br. Both the size of the halogen atoms and the reduced symmetry prevented 
us from performing core- correlation calculations for nonidentity SAr2 reactions, except for 
the F/Cl nonidentity case. 

From Table 2, it can be seen that the final Wl' and W2h energy values for the iden- 
tity reactions are very close to each other. Considering the very close agreement between 
Wl' and W2h results, the conclusion is warranted that the results from Wl' theory can 
be used as reference values to compare the results from other methods when W2h results 
are not available. As a general observation, the core contributions for the transition state 
structures are noticeably larger than for the ion-molecule complexes (see Table 2). Although 
the core-correlation contribution is small in absolute terms, its relative contribution to the 
likewise small overall barrier heights can be substantial. For example, it is 0.36 kcal/mol 
for [F- ■ -CHs" ■ -F]^, while the total energy is -0.37 kcal/mol. Likewise, the core correlation 
contributions for [CI- ■ -CHs- ■ -CI]" and [Br- ■ -CHs- ■ -Br]" are nearly 10% and 25% of the 
overall barrier (relative to reactants), respectively. At the W2h level, the core correlation 
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contribution to the total energy increases slightly. Scalar relativistic effects exhibit simi- 
lar trends as those of core correlation, but the effects are fairly small. Only Cl~- ■ -CHsBr, 
[CI- • -CHs" ■ -Br]^ and Br^- ■ -CHsCl exhibit noticeable scalar relativistic contributions, due 
to the presence of the heavy halogen Br. 

Among the identity reactions only [F- ■ -CHs- ■ -F]" has a transition state below the reac- 
tants energy level. In the nonidentity case all the transition state structures lie below the 
react ants. 

The computed final heats of formation (AHJ) of CH3X in kcal/mol are compared with 
experiment in Table 3. Both Wl' and W2h values are presented after accounting for ZPVEs 
and thermal corrections calculated at the B3LYP/cc-pVTZ+l level. At this level the ZPVEs, 
after scaling by 0.985jii are found to be 24.16, 23.26 and 22.89 kcal/mol, respectively, for 
CH3F, CH3CI and CHsBr. The corresponding thermal corrections are -1.92, -1.89 and -3.68 
kcal/mol. The computed AHJ value for CH3CI lies within the experimental error bar: the 
experimental value for CH3F is a crude estimate (±7 kcal/mol) and our computed value is 
certainly more reliable. Our calculated value for CH3Br is slightly outside the experimental 
error bar: some of the discrepancy could be due to the limitations of the scalar relativistic 
treatment. As shown by Bauschlicher,^ the simple DMV correction starts to exhibit minor 
deficiencies for third-row compounds; for first- and second-row compounds, it is in excellent 
agreement with more rigorous treatments.i2liil 

Also included in Table 3 are calculated electron affinities of X (X = F, CI and Br) in eV 
together with experimental results. Using a similar approach, but with even larger spdfghi 
basis sets as well as full CI corrections, we were ablell to reproduce the experimental EAs 
of the 1st- and 2nd-row atoms to within ±0.001 eV on average. The presently calculated 
W2h results of F and CI EAs differ by only about 0.001 eV from these benchmark values 
(EA(F)=3.403 eV and EA(C1)=3.611 eV), and the W2h results for F, CI and Br are all 
within 0.003 eV of experiment. Although the Wl' values differ about 0.01 eV for F and CI 
and 0.02 eV for Br, this is comparable to the W1/W2 target accuracy (0.25 kcal/mol on 
average). The performance of the Wl' and W2h methods for the reactants and products is 
obviously encouraging for the study of the problem at hand. 
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C. Siv2 Reactions 



The reaction mechanism with the double-well potential energy surface for the gas-phase 
Sn2 reactions is shown in Figure 1. Obviously, the energy profile is symmetric for the 
identity reactions (Figure la), and asymmetric for the nonidentity reactions (Figure lb). 
The complexation energy (AHcomp), central barrier (Alficent), ^-nd overall activation barrier 
relative to the separated reactants {AH^ovr) are defined in Figure 1. In the nonidentity case, 
the following additional quantities are defined in Figure lb: overall enthalpy change for the 
reaction (AHgyr) and the central enthalpy difference AHcent between product and reactant 
ion-molecule complexes, 3 and 1. 

D. Identity Reactions 

Complexation energies {AHcomp)-i overall barrier heights {Alt'ovr) and central barriers 
{Ati^cent) obtained from Wl' and W2h methods are compared in Table 4 with DFT, Gn, 
and CBS-QB3 methods together with available experimental values. 

It should be emphasized that the experimental data for the SAr2 reactions are insufficient 
and the available data are subject to large uncertainties. Therefore, it would be appropriate 
to analyze the performance of various methods with respect to Wn methods. First of all, 
note that the mPWlK/6-31+G* AHcomp (13.55 kcal/mol) for the F/F case is very close to 
the Wl' and W2h results (13.66 and 13.72 kcal/mol, respectively). mPWlK/cc-pVDZ(+X) 
and mPWlK/cc-pVTZ(-|-X) methods however predict lower AHcomp values. In fact, the 
B3LYP, B97, HCTH-120, mPWlPWQl, mPWH&HPW91 and mPWlK methods all predict 
roughly 1 kcal/mol lower complexation energies, while BH&HLYP and B97-1 agree well with 
Wl' and W2h. (The HCTH AHcomp is much lower than the others, vide infra.) G2 and CBS- 
QB3 values are close to the Wn results while the G3 method predicts higher complexation 
energy compared to the Wn methods. Inclusion of diffuse functions for the Gn and CBS-QB3 
reference geometries (i.e. Gn(-|-) and CBS-QB3(-|-)) increases the AHcomp value by 0.4-0.6 
kcal/mol. 

A comparison of overall barrier heights is presented in the third column of Table 4. Both 



12 



Wl' (-0.37 kcal/mol) and W2h (-0.34 kcal/mol) theories predict negative barrier heights in 
the F/F case and the values are very close. The CBS-QB3(-|-) result is in excellent agreement 
therewith; all Gn theories predict barrier heights that are lower by 1 kcal/mol, with further 
lowering seen at the G2(-|-) and G3(+) levels. Among the DFT methods considered, only 
mPWlK/6-3H-G* and HCTH/cc-pVDZ(-|-X) fortuitously predict overall barrier heights 
close to the Wn results: basis set extension for mPWlK leads to positive overall barrier 
heights, which are likewise found for the "half and half functionals. B3LYP, B97(-l) and 
HCTH-120 all significantly underestimate the barrier, mPWlPWQl to a lesser extent. 

For the chlorine identity gas-phase SAr2 reactions, fairly accurate experimental values 
are available and are presented in Table 4. The experimental values reported by Li and 
coworker^ll correspond to the standard state. Hence, thermal corrections and ZPVEs are 
subtracted from experimental values in order to compare with the "bottom of the well" 
calculated values. It is noteworthy that the Wl' (10.54 kcal/mol) and W2h (10.94 kcal/mol) 
complexation energies are in good agreement with the experimental value (10.53 kcal/mol). 
CBS-QB3 results are also in agreement with the Wn and experimental values, while those 
from DFT calculations are less satisfactory as they are about 1 kcal/mol lower. Also note 
that the Gl and G2 methods reproduce the complexation energy well, while G3 results are 
0.5 kcal/mol higher than the Wn and experimental values. 

The overall barrier height for the Cl/Cl reaction is found to be 3.07 and 2.67 kcal/mol 
at the Wl' and W2h levels of theory. Note first that the experimental value (2.90 kcal/mol) 
is very close and lies between the Wl' and W2h values. The Alfiovr value calculated at the 
mPWlK/6-31+G*, and CBS-QB3 levels of theory as well as the G2(+) value by Radom et 
a/.@ and the CCSD (T) / spdf g value by Botschwinail agree well with the W2h result. The 
mPWlK exchange-correlation functional with the cc-pVDZ(+X) and cc-pVTZ(+X) basis 
sets predict somewhat higher Alfiovr values, while the Gl, G2MP2 and G3(+) values are 
about 1 kcal/mol lower. B3LYP, B97, B97-1, and HCTH-120 all predict a negative overall 
barrier for the Cl/Cl system, in disagreement with all other methods considered and with 
experiment. BH&HLYP performs moderately well, while mPWH&HPW91 predicts a larger 
barrier height (4.50 kcal/mol) than Wn. The central barrier values presented in the last 
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column of the Table 4 reveal that the agreement between Wn theories (13.61 kcal/mol) 
and experiment (13.66 kcal/mol) is excellent. The G2MP2, G2(+) and CBS-QB3 methods 
also reproduce the central barriers very well. As expected from the overall barrier heights, 
the DFT results are less satisfactory, except for mPWlK/cc-pVDZ(+X) and mPWlK/cc- 
pVTZ(+X) which are in good agreement with the Wn values. 

For the bromine identity SAr2 reaction, Wl' theory predicts 10.03 kcal/mol for AHcomp- 
Note that the G2(+) value is in close agreement with Wl' theory. The reported experimental 
value (11.34±0.4 kcal/mol) agrees fairly well. Most DFT levels of theory considered suggest a 
complexation energy about 1 kcal/mol lower than the Wl' value, except mPWlK/6-3H-G** 
which is higher (12.78 kcal/mol, probably an artifact of the small basis set); B97 and B97- 
1 which closely bracket the Wl' value; HCTH-120 which is close to the Wl' value (see 
below); and HCTH which is 2.5 kcal/mol lower than the latter. The complexation energies 
for X~---CH3X are found to decrease in the order F > CI > Br. This trend was noted 
previously by Radom and coworkers,^! who attributed it to the electronegativities of the 
halogens. 

The overall barrier height for the Br/Br reaction is found to be 1.02 and 0.77 kcal/mol at 
the Wl' and W2h levels of theory. Of the various exchange-correlation functionals consid- 
ered, only BH&HLYP, mPWH&HPW91, and mPWlK find positive barriers (as do the Gn 
theories). It should be pointed out that the DFT results for this system display appreciable 
basis set sensitivity: for instance, the mPWlK/6-31+G* overall barrier has the wrong sign. 
It is interesting to note that the complexation energy derived from the experimental overall 
(1.73 kcal/moli) and central (11.68 kcal/molS) barrier heights is 9.95 kcal/mol while the 
reported experimental complexation energy (11.34 kcal/mol)§ is inconsistent with the de- 
rived value. In fact, the derived value is in excellent agreement with the Wl' value (10.03 
kcal/mol). This clearly suggests that the experimental data should be re-examined. 

The performance of both B97 and B97-1 for the identity reactions is quite similar to that 
of B3LYP. While the 'pure DFT' HCTH functional appears to yield markedly better over- 
all barrier heights, this comes at the expense of significantly underestimated complexation 
energies (and severely overestimated ion- molecule distances, see Supplementary Material). 
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It was previously notedE^ that HCTH severely underestimates interaction energies of H- 
bonded complexes; this was ascribed to the absence of anions and H-bonded dimers in the 
original HCTH parametrization set. A reparametrization@ against an enlarged sample of 
high-quality ab initio energies, denoted HCTH-120, eliminates this particular problem.S For 
the identity SAr2 reactions, we find that complexation energies (and ion-molecule distances) 
are dramatically improved compared to HCTH: no corresponding improvement is however 
seen for the central barrier heights, and the overall barrier heights deteriorate accordingly. 

Overall, the DFT methods are less satisfactory for barrier height calculations. Although 
the performance of mPWlK/6-3H-G* method for F/F and Cl/Cl reactions was excellent, 
it is not the ultimate low cost method for barrier heights as it has predicted a negative 
barrier for the Br/Br system. This behavior illustrates the inadequacy of the 6-31+G* basis 
set for Br: the more extended correlation consistent basis sets with the mPWlK exchange- 
correlation functional do predict the sign correctly. In addition, Gn(+) and CBS-QB3(+) 
provide an acceptable account of reaction energetics. 

E. Nonidentity Reactions 

A comparison of computed and experimental complexation energies for the nonidentity 
Si\f2 reactions is provided in Table 5. For the F~- ■ -CHsCl ion-molecule complex we could 
find a stationary point neither at the MP2/6-31G* level of theory used for the G2 and G3 
reference geometries nor at the B3LYP/6-311G(2(i, d,p) level used for the CBS-QB3 reference 
geometries; at these levels of theory, the optimization leads to Cl~- ■ -CHsF even if the initial 
geometries were chosen to correspond to F~- ■ -CHsCl. Addition of diffuse functions to the 
basis set for the reference geometry remedies the problem. Similarly, in the F/Br case only 
the Br^ • ■ -CHsF complex is found as a stationary point at the MP2/6-31G* level of theory, 
and the transition state and second ion-molecule complex only appear when diffuse functions 
are added to the basis set. Furthermore, and regardless of the basis set employed, none of 
the DFT functional except mPWlK, mPWH&HPW91, and BH&HLYP find a transition 
state or a F^ • • -CHaBr complex. (The CBS-QB3 method is not defined for Br and hence 
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no CBS-QB3 data are presented for the F/Br and CI/Br ion-molecule complexes.) Table 5 
also presents large-scale CCSD(T) energ etics for the F/Cli and CI/bH cases reported by 
Botschwina and coworkers. Available experimental values are presented at the end of the 
Table with uncertainties in parentheses. 

Examination of Table 5 indicates that the AHcomp values strongly depend on the nucle- 
ophile (Y~), decreasing in the order F~ > Cl^ > Br^. They also depend on the leaving 
group (X~), in the order CH3F < CH3CI < CHsBr. Similar observations were made earlier 
by Radom and coworkers.0 

Comparison of the complexation energies obtained from various methods with Wl' theory 
indicates that all DFT results for Cl~- • -CHsF are lower by 1 kcal/mol. The only available 
experimental valueil for Cl^- ■ -CHaF {AH° = 11.41 kcal/mol) has an uncertainty of 2.01 
kcal/mol. Comparison of this value with the calculated values suggest that more accurate 
measurements are in order. For Cl~- ■ -CHsBr and Br~- ■ -CHsCl, rather more accurate high- 
pressure mass spectrometry data are available (12.54 and 11.01 kcal/mol). The Wl' values 
(11.91 and 10.32 kcal/mol) are very close to the experimental results, considering the exper- 
imental uncertainty of 0.4 kcal/mol. While the mPWlK/6-31+G* values for Cl~- ■ -CHsBr 
and Br^- ■ -CHaCl are fortuitously within the experimental error bars, the other DFT meth- 
ods predict lower values. Also note that G2(+) predicts complexation energies close to Wl' 
and experiment for Cl^- ■ -CHaBr, while the Br~- ■ -CHsCl value is small. A complete as- 
sessment of CBS-QB3 is not possible as it could not be applied to the bromine- containing 
systems. 

Like for the identity reactions, complexation energies are significantly underestimated 
(and ion-molecule distances overestimated by up to 0.3 A: see Supplementary Material) by 
HCTH, and this problem is mostly remedied by HCTH-120. B97 and especially B97-1 appear 
to represent an improvement over B3LYP for the complexation energies. 

Calculated overall reaction enthalpies, central enthalpy differences between reactant and 
product ion-molecule complexes, overall barrier heights, and central barrier heights for the 
nonidentity SAr2 reactions are presented with available experimental results in Table 6. It 
needs to be reemphasized that all values are "bottom-of-the-well" (i.e, zero-point exclusive): 
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that is, the experimental values are presented after subtracting the ZPVEs (scaled by 0.985) 
and thermal corrections obtained using the B3LYP method. 

The overall reaction enthalpies of the three nonidentity reactions, viz F/Cl, F/Br, and 
CI/Br, calculated at the Wl' level are -32.65, -41.43 and -8.56 kcal/mol, respectively. The 
corresponding experimental values are available and are presented in Table 6. The exper- 
imental value for the F/Cl reaction (-33.34 kcal/mol) is in close agreement with the Wl' 
value. B3LYP, B97(-l), mPWlPWQl, Gn, and CCSD(T)/spdfg results are all in close 
agreement with the Wl' value, but mPWlK, "half and half and CBS-QB3 theories pre- 
dict 3-5 kcal/mol higher exothermicity. For the F/Br reaction, the mPWlK/6-3H-G* and 
mPWlPW91/cc-pVTZ(-|-X) methods yield overall reaction enthalpies which are quite close 
to the Wl' result. The experimental result (-40.20 kcal/mol) is in good agreement with the 
best calculated values considering the uncertainty of 1 kcal/mol. The Gn theories predict 
exothermicities below, and mPWlK/cc-pVnZ(-|-X) above, the Wl' value. Concerning the 
CI/Br nonidentity reaction, the reported experimental value (-6.86 kcal/mol) differs from 
the Wl' value by 2 kcal/mol. As expected, the very recent CCSD (T) / spdf gh results of 
Botschwina and coworkers^ are in close agreement with our predictions. Our results suggest 
that the CI/Br experimental data may need to be reconsidered. Note the significant basis set 
dependence in the mPWlK results, which illustrates the inadequacy of the 6-31-l-G* basis 
set. 

As a general observation, 50:50 admixture of HF exchange in the DFT theories increases 
the magnitude of the overall reaction enthalpy, and the increase is greater in BH&HLYP than 
in mPWH&HPWQl. Performance of B97 and B97-1 for the overall reaction enthalpies is 
similar to B3LYP, while HCTH and HCTH-120 represent underestimates in absolute value. 

At the Wl' level of theory, the calculated central barrier (Alficent) for the F/Cl system 
is 2.89 kcal/mol. Gn{+), CBS-QB3(+), BH&HLYP, mPWH&HPW91, and mPWlK all 
reproduce the Wl' value moderately well, while B3LYP, mPWlPW91, B97, B97-1, HCTH, 
and HCTH-120 all underestimate the central barriers. The experimental central barrier for 
the F/Cl system is nearly 4 kcal/mol higher than the calculated values. Judging from the 
performance of the various methods for the identity SAr2 reactions, it is almost certain that 
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the experimental F/Cl central barrier is in error and unambiguous new measurements are 
in order. For the F/Cl and CI/Br systems, the B3LYP, mPWlPW91, B97, B97-1, HCTH, 
and HCTH-120 central barriers are all underestimated, while these exchange-correlation 
functionals find no barrier at all for the F/Br case. Like for the identity case, mPWlK/cc- 
pVTZ(+X) and BH&HLYP/cc-pVTZ(+X) central barriers agree well with the benchmark 
ab initio values, although the basis set sensitivity of particularly the CI/Br results argues 
against using small basis sets like 6-31+G*. 

Several studies have reported experimental overall barrier heights (Alfiovr), but only for 
the CI/Br system, and the experimental data range from -0.61 to -1.83 kcal/mol. To our 
knowledge, no experimental data are available for the F/Cl and F/Br systems. For the CI/Br 
system, the theoretical values span a range from -1.17 to -6.60 kcal/mol. Nevertheless, it 
is worth noting that the Wl' value (-1.82 kcal/mol) for the CI/Br system is in excellent 
agreement with the experimental overall barrier height reported by Caldwell and coworkers!! 
(-1.83 kcal/mol, after accounting for ZPVE and thermal corrections). Some caution should 
be exercised as the Wl' value does not include the core correlation contribution. Also note 
that G2 theory (-1.82 kcal/mol) reproduces the Wl' value very well. Except for mPWlK/cc- 
pVTZ(+X), mPWH&HPW91/cc-pVTZ(+X), and BH&HLYP/cc-pVTZ(+X), all the DFT 
methods perform poorly, consistent with the preceding discussion. 



IV. CONCLUSIONS 

A benchmark study using the Wl' and W2h methods has been carried out for the po- 
tential surface of the gas-phase 3^2 reactions -|- CH3X — > CH3Y -|- X^. A number 
of more approximate (and less expensive) methods — both compound models (like G2/G3 
theory and CBS-QB3) and density functional methods — have been applied in an attempt 
to assess their performance for barrier heights in Sn^ reactions. We arrive at the following 
conclusions. 

(1) Our best calculations are in excellent agreement with experiment for the AH J values 
of the methyl halides (where available) and the electron affinities of the halogens. Where 
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accurate experimental data are available for the title reactions (e.g. for the Cl/Cl case), our 
best calculations agree with experiment to within overlapping uncertainties. Our calculations 
however suggest that more reliable experimental data are in order for most of the reactions 
considered. 

(2) The nonidentity S]sf2 reactions and F/F identity reaction possess transition state struc- 
tures below the reactants energy while Cl/Cl and Br/Br transition structures are above the 
reactants energy. The complexation energies for identity Sjv2 reactions are found to increase 
in the order Br < CI < F while the barrier heights follow the order F ^ Br < CI. The 
complexation energies for the nonidentity reactions indicate that the AHcomp strongly 
depend on the nucleophile and leaving group. 

(3) The B3LYP, and to a lesser extent, mPWlPW91 exchange-correlation functionals sys- 
tematically underestimate barrier heights and, in the F/Br case, are not even able to locate 
the correct stationary points on the potential surface. The latter problem is remedied by 
using the corresponding 'half-and-half functionals BH&HLYP and mPWH&HPW91, which 
however appear to somewhat overcorrect the barrier height. The B97 and B97-1 function- 
als perform similarly to B3LYP for the problem under study. The 'pure DFT' HCTH and 
HCTH-120 functionals both underestimate central barrier heights; HCTH in addition under- 
estimates complexation energies (and severely overestimates ion-molecule distances), which 
are however well reproduced by HCTH-120. Overall, the mPWlK functional appears to 
put in the best performance of all DFT methods considered, especially when using extended 
basis sets. 

(4) The performance of G2(-|-), G3(-|-), and CBS-QB3(-|-) methods for the energetics still 
appears to be superior to the DFT methods. (The '(-l-)' stands for the addition of diffuse 
functions to the basis set used in obtaining the reference geometries; this is mandatory to 
get transition states at all in the F/Br and CI/Br cases.) The limitations for transition 
states of the B3LYP exchange-correlation functional suggest its replacement — at least for 
kinetics applications — by mPWlK in thermochemistry methods that employ DFT reference 
geometries, e.g. G3B3, CBS-QB3, and Wl theory. (In addition, larger basis sets than 6- 
31-l-G* should definitely be considered for the Br compounds.) 
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The present calculations illustrate the power of state-of-the-science theoretical methods 
in providing both qualitative and quantitative information regarding the reaction energetics. 
In the absence of accurate experimental data, our high quality results should be useful to 
future experimental and theoretical studies. 
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TABLES 

TABLE 1. Components of computed electron affinities of X and total atomization energies (kcal 
mol-i) of CH3X (X = F, CI and Br). 
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F- 


30.08 


44.71 


4.15 


0.17 


-0.39 


-0.26 


78.46 


ci- 


58.33 


23.82 


2.26 


0.02 


-0.84 


-0.34 


83.25 


Br- 


58.31 


21.49 


1.91 


0.31 


-3.51 


-0.90 


77.62 


CH3F 


319.82 


96.89 


5.34 


1.14 


-0.47 


-0.37 


422.34 


CH3CI 


303.90 


86.03 


5.28 


1.21 


-0.93 


-0.42 


395.07 


CH:5Br 


293.11 


85.15 


5.21 


1.19 


-3.60 


-0.79 


.380.91 
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TABLE 2. Components of relative energies 


(kcal mol 


of ion- 


■molecule complexes 


and transi- 


tion state structures with respect to reactants. 












SCF 


CCSD 


(T) 


Core 


Scalar Rel. 


Final 


Species 


limi( 


limit 


limit 


corr. 






Wl' 














• CHaF 


-10.63 


-2.54 


-0.56 


0.09 


-0.02 


-13.66 


[F-CH3--F]- 


8.23 


-6.23 


-2.65 


0.36 


-0.08 


-0.37 


Cr- • -CHaCl 


-8.08 


-1.93 


-0.62 


0.08 


0.01 


-10.54 


[CI- • -CHa- • -CI] 


7.81 


-2.62 


-2.37 


0.35 


-0.10 


3.07 


Br~- ■ CHaBr 


-7.55 


-1.92 


-0.65 




0.09 


-10.03" 


[Br- ■ -CHa- • -Br] 


5.71 


-2.43 


-2.36 


0.26 


-0.16 


1.02 


W2h 














F--CHaF 


-11.16 


-2.17 


-0.45 


0.07 


-0.02 


-13.72 


Pi — 1 ✓"NTT 1 — lT — 

[F- • -CHa- • -F] 


8.58 


-6.61 


-2.58 


0.35 


-0.08 


-0.34 


CI • • -CHaCl 


-8.01 


-2.39 


-0.65 


0.10 


0.01 


-10.94 


[CI- • -CHa- ■ -CI] 


8.55 


-3.63 


-2.46 


0.33 


-0.11 


2.67 


r -f — V y^TT "T — \ 1 

[Br- ■ -CHa- • Br] 


6.57 


-3.41 


-2.38 


0.21 


-0.21 


0.77 


Wl' 














F"- • CHaCl 


-13.67 


-1.16 


-0.76 


0.18 


-0.02 


-15.43 


[F- • -CHa- • -CI]- 


-10.39 


-0.84 


-1.60 


0.36 


-0.07 


-12.54 


Cr- • -CHaF 


-7.08 


-2.03 


-0.42 


0.03 


0.00 


-9.51 


F~- • CHaBr 


-15.08 


-1.24 


-0.72 




0.03 


-17.01" 


Ft? r^TT t3t-1 — 
[r • • -Uxla- • -JDrJ 


-io.oo 
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-1.43 
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Br-- • -CHaF 


-6.63 


-1.69 


-0.24 




0.05 


-8.51" 


Cr - CHaBr 


-8.89 


-2.02 


-0.59 




-0.41 


-11.91" 


[CI- - -CHa- - -Br]- 


2.24 


-2.69 


-2.28 




-0.57 


-3.30" 


Br----CH3C1 


-7.38 


-2.03 


-0.45 




-0.46 


-10.32" 


" Wl' — Core Correlation 
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TABLE 3. Calculated and experimental heats of formation (kcal mol ^) of CH3X and electron 
affinities (eV) of X (X = F, CI and Br) 

Heats of Formation Electron Affinity 



Species 


Wl' 


W21i 


Experiment 


Species 


Wl' 


W2h 


Experiment 


CH3F 


-57.06 


-56.75 


-56(7)i 


F 


3.411 


3.402 


3.401 190(4)S 


CH3CI 


-20.14 


-20.23 


-20.00(50)11 


CI 


3.627 


3.610 


3.612 69(6)ii 


CHsBr 


-8.50 


-8.63 


-8.20(19)§ 


Br 


3.383 


3.366 


3.363 583(44)11 
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TABLE 4. Comparison of complexation energies (A^fcomp) of the ion-molecule complexes, over- 
all barrier heights relative to reactants (Alfiovr), and central barriers (Alficent) of identity SAr2 
reactions, X~ + CH3X — > XCH3 + X~, calculated at various levels of theory. All values in 



kcal/mol. 



X 


Method 








F 


Wl' 


13.66 


-0.37 


13.29 




W2h 


13.72 


-0.34 


13.38 




B3LYP/cc-pVTZ{+X) 


12.72 


-2.58 


10.15 




BH&HLYP/cc-pVTZ(+X) 


13.22 


1.31 


14.53 




mPWlPW91/cc-pVTZ(+X) 


12.49 


-0.95 


11.55 




mPWH&HPW91/cc-pVTZ(+X) 


12.77 


2.60 


15.38 




B97/cc-pVDZ(+X) 


12.48 


-2.47 


10.01 




B97-l/cc-pVDZ(+X) 


13.21 


-3.29 


9.92 




HCTH/cc-pVDZ(+X) 


9.87 


-0.60 


9.27 




HCTH-120/cc-pVDZ(+X) 


12.39 


-4.20 


8.18 




mPWlK/6-31+G* 


13.55 


-0.30 


13.26 




mPWlK/cc-pVDZ{+X) 


12.63 


0.36 


13.00 




mPWlK/cc-pVTZ(+X) 


12.66 


1.69 


14.36 




Gl 


13.01 


-1.37 


11.64 




G2 


13.34 


-1.15 


12.19 




G2MP2 


13.41 


-0.63 


12.78 




G3 


14.23 


-1.97 


12.26 




CBS-QB3 


13.46 


-0.85 


12.61 




G3(+) 


14.59 


-2.68 


11.90 




CBS-QB3(+) 


14.15 


-0.52 


13.63 




G2(+)'' 


13.81 


-1.86 


11.95 




CCSD{T) /spdf 


13.73 


-0.92 


12.81 


CI 


Wl' 


10.54 


3.07 


13.61 




W2h 


10.94 


2.67 


13.61 




B3LYP/cc-pVTZ{+X) 


9.50 


-0.48 


9.02 




BH&HLYP/cc-pVTZ(+X) 


9.67 


3.17 


12.84 




mPWlPW91/cc-pVTZ(+X) 


9.59 


1.23 


10.82 




mPWH&HPW91/cc-pVTZ(+X) 


9.69 


4.50 


14.19 




B97/cc-pVDZ(+X) 


10.10 


-0.66 


9.44 




B97-l/cc-pVDZ(+X) 


10.74 


-1.46 


9.28 




HCTH/cc-pVDZ(+X) 


7.91 


1.45 


9.36 




HCTH-120/cc-pVDZ(+X) 


9.96 


-1.93 


8.03 




mPWlK/6-31+G* 


9.75 


3.20 


12.95 




mPWlK/cc-pVDZ{+X) 


9.65 


3.63 


13.28 




mPWlK /cc-pVTZ(+X) 


9.64 


3.66 


13.30 




Gl 


10.52 


1.79 


12.31 




G2 


10.77 


3.06 


13.83 
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G2MP2 


10.89 


2.74 


13.63 


G3 


11.15 


1.79 


12.95 


CBS-QB3 


10.65 


2.47 


13.12 


G3(+) 


11.04 


1.80 


12.83 


CBS-QB3(+) 


10.69 


2.40 


13.09 


G2(+)'' 


10.71 


3.01 


13.72 


CCSB{T) /spdfgl' 




2.75 




Experiment 


10.53(40)'* 


2.90"= 


13.66(2.01; 


Wl' — core 


10.03 


1.02' 


10.79 


W2h 




0.77 




B3LYP/cc-pVTZ(+X) 


9.06 


-2.42 


6.64 


BH&HLYP/cc-pVTZ{+X) 


9.04 


1.25 


10.29 


mPWlPW91 /cc-pVTZ(+X) 


9.21 


-1.03 


8.18 


mPWH&HPW91 /cc-pVTZ(+X) 


9.19 


2.22 


11.40 


B97/cc-pVDZ(+X) 


9.62 


-2.29 


7.33 


B97-l/cc-pVDZ(+X) 


10.24 


-3.02 


7.22 


HCTH/cc-pVDZ(+X) 


7.56 


-0.70 


6.86 


HCTH-120/cc-pVDZ(+X) 


9.73 


-4.06 


5.68 


mPWlK/6-31+G* 


12.78 


-1.95 


10.83 


mPWlK/cc-pVDZ(+X) 


9.34 


0.68 


10.02 


mPWlK/cc-pVTZ(+X) 


9.16 


1.38 


10.54 


Gl 


9.68 


1.11 


10.78 


G2 


9.85 


1.52 


11.38 


G2MP2 


9.83 


1.83 


11.66 


G2(+)" - ECP 


10.17 


1.48 


11.65 


Experiment 


11.34(40)'' 


1.73s 


11.68' 



" G2(+) values are from Ref.@; ^ CCSD{T) /spdfg values are from Ref.i; " ReiE 
Experimental values: '^From Ref.0, ^ From Ref.S, ^ From Ref.0, 9 From Ref.ii, ^ From 
Ref.i 

* Core contribution included 
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TABLE 5. Comparison of complexation energies {AHcomp, kcal/mol) of the ion-molecule 
complexes for the nonidentity SAr2 reactions, calculated at various levels of theory. 



Method 


F"- • -CHaCl 


CI"- • CHsF 


F~- ■ CHsBr 


Br~- • -CHsF 


CI" - • CHsBr 


Br"- ■ CHsCl 


Wl' — core 


15.43" 


9.51"^ 


17.01 


8.51 


11.91 


10.32 


B3LYP/cc-pVTZ(+X) 


15.37 


8.09 




7.18 


10.24 


8.42 


BH&HLYP/cc-pVTZ(+X) 


15.39 


8.42 


16.65 


7.45 


10.23 


8.56 


mPWlPWQl /cc-pVTZ(+X) 


15.06 


8.13 




7.27 


10.32 


8.54 


mPWH&HPW91 /cc-pVTZ(+X) 


15.03 


8.32 


16.28 


7.42 


10.30 


8.65 


B97/cc-pVDZ(+X) 


15.43 


8.49 




7.69 


10.66 


9.13 


B97-l/cc-pVDZ(+X) 


16.26 


9.05 




8.21 


11.32 


9.74 


HCTH/cc-pVDZ{+X) 


12.66 


6.56 




5.82 


8.56 


7.00 


HCTH-120/cc-pVDZ(+X) 


15.42 


8.48 




7.69 


10.78 


9.01 


mPWlK/6-31+G* 


15.32 


8.79 


17.63 


9.12 


12.97 


10.56 


mPWlK/cc-pVDZ (+X) 


14.52 


8.46 


16.25 


7.62 


10.41 


8.68 


mPWlK/cc-pVTZ{+X) 


14.97 


8.25 


16.30 


7.37 


10.27 


8.60 


Gl 




9.58 




8.35 


11.17 


9.10 


G2 




9.68 




8.42 


11.35 


9.36 


G2MP2 




9.71 




8.40 


11.37 


9.45 


G3 




10.03 










CBS-QB3 




9.33 










G3(+) 


16.34 


9.97 










CBS-QB3{+) 


15.85 


9.51 










G2(+)'' 


15.62 


9.64 


16.74 


8.56 


11.47 


9.64 


CCSD(T)/large= 


16.07 


9.75 






11.31 


9.71 


Experiment. 




11.41(2.01)'* 






12.54(40)*= 


11.01 (40)*= 



Core contribution included 



G2(-|-) values are from Ref.0, 
" F/Cl: CCSD {T)/spdfg values from Ref.il; CI/Br: CCSD{T) /spdfgh values from Ref.S 
Experimental values: '^From Ref.S, ^ From Ref.0 
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TABLE 6. Comparison of overall reaction enthalpies {^Hovr)-, central enthalpy differences 
between reactant and product ion-molecule complexes (Ai/cent)) overall barrier heights {Alfiover) 
and central barrier heights (Ai/^cent) for exothermic + CH3X — > YCH3 + X~ reactions, 



calculated at various levels of theory. All values in kcal/mol. 



Y,X 


Method 






over 




F/Cl 


Wl' 


-32.65 


-26.73 


-12.54 


2.89 




B3LYP/cc-pVTZ(+X) 


-32.77 


-25.49 


-14.69 


0.67 




BH&HLYP/cc-pVTZ(+X) 


-37.02 


-30.05 


-12.86 


2.53 




mPWlPW91/cc-pVTZ(+X) 


-33.08 


-26.15 


-13.43 


1.63 




mPWH&HPW91/cc-pVTZ(+X) 


-36.45 


-29.73 


-11.51 


3.52 




B97/cc-pVDZ(+X) 


-32.90 


-25.95 


-14.70 


0.73 




B97-l/cc-pVDZ(+X) 


-33.13 


-25.92 


-15.60 


0.66 




HCTH/cc-pVDZ(+X) 


-30.77 


-24.67 


-11.95 


0.71 




HCTH-120/cc-pVDZ(+X) 


-30.58 


-23.64 


-15.14 


0.27 




mPWlK/6-31+G* 


-36.59 


-30.07 


-13.02 


2.30 




mPWlK/cc-pVDZ (+X) 


-34.74 


-28.68 


-11.95 


2.57 




mPWlK/cc-pVTZ(+X) 


-35.50 


-28.78 


-11.97 


3.01 




Gl 


-30.62 










G2 


-31.59 










G2MP2 


-32.25 










G3 


-33.00 










CBS-QB3 


-35.21 










G3(+) 


-32.86 


-26.50 


-14.04 


2.30 




CBS-QB3(+) 


-35.15 


-28.81 


-13.77 


2.07 




G2(+)'' 


-31.44 


-25.46 


-12.63 


2.98 




CCSD{T)/spdfg>' 


-32.34 


-26.36 


-11.84 


3.89 




Experiment 


-33.34(72)= 






7.52(1.20)'* 


F/Br 


Wl' — core 

B3LYP /cc-pVTZ(+X) 


-41.43 

-40.78 


-32.93 


-16.37 


0.64 




BH&HLYP/cc-pVTZ(+X) 


-46.13 


-36.93 


-16.33 


0.32 




mPWlPW91 /cc-pVTZ (+X) 


-41.65 










mPWH&HPW91/cc-pVTZ(+X) 


-45.90 


-37.05 


-15.42 


0.85 




B97/cc-pVDZ(+X) 


-40.41 










B97-l/cc-pVDZ(+X) 


-40.68 










HCTH/cc-pVDZ(+X) 


-38.18 










HCTH-120/cc-pVDZ(+X) 


-37.96 










mPWlK/6-31+G* 


-42.35 


-33.84 


-16.92 


0.71 




mPWlK/cc-pVDZ (+X) 


-44.76 


-36.13 


-16.08 


0.18 




mPWlK/cc-pVTZ(+X) 


-44.70 


-35.77 


-15.76 


0.54 




Gl 


-40.87 










G2 


-40.02 
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CI/Br 



G2MP2 


-40.35 








G2(+)'^ 


-39.47 


-31.29 


-15.90 


0.84 


Experiment 


-40.20(96)'= 








Wl' — core 


-8.56 


-6.97 


-1.82 


8.61 


B3LYP/cc-pVTZ(+X) 


-8.01 


-6.19 


-5.25 


4.99 


BH&HLYP/cc-pVTZ(+X) 


-9.11 


-7.43 


-2.15 


8.08 


mPWlPW91/cc-pVTZ(+X) 


-8.57 


-6.79 


-3.99 


6.33 


mPWH&HPW91/cc-pVTZ(+X) 


-9.45 


-7.79 


-1.17 


9.14 


B97/cc-pVDZ(+X) 


-7.52 


-5.99 


-5.06 


5.60 


B97-l/cc-pVDZ(+X) 


-7.54 


-5.97 


-5.83 


5.48 


HCTH/cc-pVDZ(+X) 


-7.41 


-5.85 


-3.18 


5.38 


HCTH-120/cc-pVDZ(+X) 


-7.38 


-5.62 


-6.52 


4.26 


mPWlK/6-31+G* 


-5.75 


-3.34 


-3.12 


9.85 


mPWlK/cc-pVDZ(+X) 


-10.02 


-8.30 


-6.60 


3.81 


mPWlK/cc-pVTZ(+X) 


-9.20 


-7.53 


-1.88 


8.38 


Gl 


-10.25 


-8.18 


-3.45 


7.72 


G2 


-8.43 


-6.44 


-1.82 


9.53 


G2MP2 


-8.11 


-6.18 


-1.67 


9.70 


G2(+)<' 


-8.04 


-6.21 


-1.71 


9.76 


CCSD {T)/spdfgh'' 


-8.53 


-6.93 


-2.33 


8.98 


Experiment 


-6.86(72)'= 




-1.83(5)"= 





-1.69(33)-^ 
-1.529 
-Lll*" 
-0.61' 

" G2(+) values are from Ref.Ei 
^ Ref.i (F/Cl) and Ref.il (CI/Br) 
^' Table 2 of Ref.i 

Experimental values: '^From Ref-S, From Ref.S, ^ From Ref.i, From Ref.ii, ^ From Ref.S, 
From RefE, ' From RefS 
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FIGURES 

FIG. 1. Schematic representation of potential energy surface for the (a) identity and (b) non- 
identity Siv2 reactions 
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TABLES 



TABLE I. Calculated and experimental geometries (A , degree) of CH3X (X=F, CI, and Br). 



Species 


Method/Basis Set 


r(G-X) 


r(C-H) 


ZXCH 


CH3F 


CCSD(T)/cc-pVQZ+l 


1.382 


1.089 


108.9 




B3LYP/cc-pVTZ(+X) 


1.392 


1.090 


108.7 




BH&HLYP/cc-pVTZ(+X) 


1.372 


1.082 


108.8 




mPWlPW91/cc-pVTZ(+X) 


1.380 


1.090 


108.9 




mPWH&HPW91/cc-pVTZ(+X) 


1.362 


1.083 


109.1 




B97/cc-pVDZ(+X) 


1.399 


1.102 


108.7 




B97-l/cc-pVDZ(+X) 


1.398 


1.102 


108.7 




HCTH/cc-pVDZ(+X) 


1.400 


1.103 


108.8 




HCTH-120/cc-pVDZ(+X) 


1.402 


1.104 


108.7 




mPWlK/6-31+G* 


1.374 


1.087 


108.7 




mPWlK/cc-pVDZ(+X) 


1.378 


1.093 


108.8 




mPWlK/cc-pVTZ(+X) 


1.367 


1.085 


109.0 




MP2/6-31G* 


1.390 


1.092 


109.1 




MP2/6-31+G*" 


1.407 


1.090 


108.0 




B3LYP/6-311G(2d,d,p) 


1.388 


1.093 


109.2 




B3LYP/6-311+G(2d,d,p) 


1.396 


1.091 


108.6 




Expt.'' 


1.383 


1.086 


108.8 


CH3CI 


CCSD(T)/cc-pVQZ+l 


1.783 


1.085 


108.4 




B3LYP/cc-pVTZ(+X) 


1.796 


1.085 


108.3 




BH&HLYP/cc-pVTZ(+X) 


1.779 


1.077 


108.4 




mPWlPW91 /cc-pVTZ(+X) 


1.776 


1.085 


108.5 




mPWH&HPW91/cc-pVTZ(+X) 


1.761 


1.079 


108.6 




B97/cc-pVDZ(+X) 


1.808 


1.097 


108.1 




B97-l/cc-pVDZ(+X) 


1.808 


1.097 


108.0 




HCTH/cc-pVDZ(+X) 


1.793 


1.098 


108.4 




HCTH-120/cc-pVDZ(+X) 


1.796 


1.098 


108.4 




mPWlK/6-31+G* 


1.772 


1.084 


108.8 




mPWlK/cc-pVDZ(+X) 


1.768 


1.090 


108.5 




mPWlK/cc-pVTZ(+X) 


1.765 


1.081 


108.6 




MP2/6-31G* 


1.777 


1.088 


108.9 




MP2/6-31+G*" 


1.780 


1.089 


108.9 




B3LYP/6-311G(2d,d,p) 


1.803 


1.087 


108.2 




B3LYP/6-311+G(2d,d,p) 


1.803 


1.087 


108.2 




Expt.'^ 


1.776 


1.085 


108.6 


CHgBr 


CCSD(T)/cc-pVQZ+l 


1.944 


1.084 


107.8 




B3LYP/cc-pVTZ(+X) 


1.957 


1.083 


107.6 




BH&HLYP/cc-pVTZ(+X) 


1.938 


1.076 


107.8 




mPWlPW91 /cc-pVTZ(+X) 


1.936 


1.084 


107.9 




mPWH&HPW91/cc-pVTZ(+X) 


1.917 


1.078 


108.0 




B97/cc-pVDZ(+X) 


1.959 


1.097 


107.6 



B97-l/cc-pVDZ(+X) 1.959 1.097 107.6 

HCTH/cc-pVDZ(+X) 1.952 1.097 107.9 

HCTH-120/cc-pVDZ(+X) 1.956 1.097 107.8 

mPWlK/6-31+G* 1.925 1.083 108.2 

mPWlK/cc-pVDZ(+X) 1.927 1.089 107.9 

mPWlK/cc-pVTZ(+X) 1.922 1.079 108.0 

MP2/6-31G* 1.947 1.086 107.9 

MP2/6-31+G*" 1.954 1.088 108.0 

Expt.'^ 1.934 1.082 107.7 

« Prom Glukhovtsev, M.N.; Pross, A.; Radom, L.; J. Am. Chem. Soc. 1995, 117, 2024. 
Experimental values: 

Prom Egawa, T.; Yamamoto, S.; Nakata, M.; Kuchitsu, K.; J. Mol. Struct.1987, 156, 213. 

Prom Jensen, T.; Brodersen, S.; Guelachvili, G.; J. Mol.Spectrosc. 1981, 88, 378. 
^ Prom Graner, G.; J. Mol.Spectrosc. 1981, 90, 394. 



TABLE II. Geometries (A , degree) of ion-molecule complexes X • • -CHaX (X=F, CI, and Br) 



of the SN2 identity reactions. 



Species 


Method/Basis Set 


r(X...C) 


r(C-X) 


r(C-H) 


ZHCX 


X=F 


CCSD(T)/cc-pVQZ+l 


2.494 


1.432 


1.082 


109.1 




B3LYP/cc-pVTZ(+X) 


2.588 


1.447 


1.082 


108.4 




BH&HLYP/cc-pVTZ(+X) 


2.581 


1.418 


1.075 


108.8 




mPWlPW91/cc-pVTZ(+X) 


2.575 


1.428 


1.083 


108.9 




mPWH&HP W9 1 /cc-pVTZ(+X) 


2.567 


1.404 


1.077 


109.2 




B97/cc-pVDZ(+X) 


2.643 


1.453 


1.095 


108.3 




B97-l/cc-pVDZ(+X) 


2.613 


1.453 


1.095 


108.3 




HCTH/cc-pVDZ(+X) 


2.794 


1.453 


1.096 


108.4 




HCTH-120/cc-pVDZ(+X) 


2.668 


1.462 


1.096 


108.2 




mPWlK/6-31+G* 


2.572 


1.421 


1.080 


108.5 




mPWlK/cc-pVDZ(+X) 


2.579 


1.425 


1.087 


108.8 




mPWlK/cc-pVTZ(+X) 


2.571 


1.411 


1.079 


109.1 




MP2/6-31G* 


2.426 


1.439 


1.083 


109.5 




MP2/6-31+G*" 


2.628 


1.456 


1.084 


107.7 




B3LYP/6-311G(2d,d,p) 


2.420 


1.453 


1.084 


110.0 




B3LYP/6-311+G(2d,d,p) 


2.565 


1.455 


1.084 


108.2 


X=C1 


CCSD(T)/cc-pVQZ+l 


3.123 


1.860 


1.080 


107.8 




B3LYP/cc-pVTZ(+X) 


3.191 


1.846 


1.080 


107.9 




BH&HLYP/cc-pVTZ(+X) 


3.203 


1.818 


1.073 


108.2 




mPWlPW91/cc-pVTZ(+X) 


3.163 


1.817 


1.081 


108.5 




mPWH&HPW91/cc-pVTZ(+X) 


3.166 


1.796 


1.075 


108.7 




B97/cc-pVDZ(+X) 


3.198 


1.855 


1.094 


107.7 




B97-l/cc-pVDZ(+X) 


3.168 


1.857 


1.093 


107.6 




HCTH/cc-pVDZ(+X) 


3.375 


1.838 


1.094 


108.1 




HCTH-120/cc-pVDZ(+X) 


3.241 


1.848 


1.095 


107.9 




mPWlK/6-31+G* 


3.181 


1.809 


1.080 


108.7 




mPWlK/cc-pVDZ(+X) 


3.182 


1.804 


1.086 


108.6 




mPWlK/cc-pVTZ(+X) 


3.167 


1.801 


1.077 


108.7 




MP2/6-31G* 


3.158 


1.812 


1.084 


109.0 




MP2/6-31+G*" 


3.270 


1.810 


1.085 


108.8 




B3LYP/6-311G(2d,d,p) 


3.146 


1.861 


1.082 


107.6 




B3LYP/6-311+G(2d,d,p) 


3.187 


1.854 


1.083 


107.8 


X=Br 


CCSD(T)/cc-pVQZ+l 


3.277 


2.023 


1.079 


106.9 




B3LYP/cc-pVTZ(+X) 


3.321 


2.016 


1.079 


106.9 




BH&HLYP/cc-pVTZ(+X) 


3.361 


1.980 


1.072 


107.4 




mPWlPW91/cc-pVTZ(+X) 


3.290 


1.982 


1.080 


107.5 




mPWH&HPW91/cc-pVTZ(+X) 


3.313 


1.953 


1.074 


107.9 




B97/cc-pVDZ{+X) 


3.320 


2.013 


1.093 


107.0 




B97-l/cc-pVDZ(+X) 


3.281 


2.014 


1.093 


106.9 




HCTH/cc-pVDZ(+X) 


3.471 


2.007 


1.094 


107.2 




HCTH-120/cc-pVDZ(+X) 


3.332 


2.022 


1.094 


106.9 



mPWlK/6-31+G* 3.215 

mPWlK/cc-pVDZ(+X) 3.304 

mPWlK/cc-pVTZ(+X) 3.312 

MP2/6-31G* 3.196 

MP2/6-31+G*" 3.395 



1.960 1.079 108.0 

1.968 1.086 107.6 

1.960 1.076 107.8 

1.992 1.083 107.7 

1.988 1.084 107.8 



« From Glukhovtsev, M.N.; Pross, A.; Radom, L.; J. Am. Chem. Soc. 1995, 117, 2024. 



TABLE III. Geometries (A , degree) of ion-molecule complexes, Y • • •CH3X, of the SN2 



non-Identity reactions (complex 1). 



Y/X 


Method/Basis Set 


r(Y...C) 


r(C-X) 


r(C-H) 


ZHCX 


F/Cl 


B3LYP/cc-pVTZ(+X) 


2.443 


1.908 


1.076 


106.1 




BH&HLYP/cc-pVTZ (+X) 


2.494 


1.853 


1.070 


107.4 




mPWlPW91/cc-pVTZ(+X) 


2.462 


1.861 


1.078 


107.3 




mPWH&HPW91/cc-pVTZ(+X) 


2.489 


1.825 


1.072 


108.1 




B97/cc-pVDZ(+X) 


2.467 


1.913 


1.090 


106.0 




B97-l/cc-pVDZ(+X) 


2.445 


1.916 


1.090 


105.9 




HCTH/cc-pVDZ(+X) 


2.565 


1.903 


1.091 


106.5 




HCTH-120/cc-pVDZ(+X) 


2.459 


1.931 


1.090 


105.6 




mPWlK/6-31+G* 


2.487 


1.846 


1.077 


107.8 




mP WIK /cc-pVDZ (+X) 


2.495 


1.837 


1.083 


107.8 




mPWlK/cc-pVTZ(+X) 


2.486 


1.833 


1.074 


108.0 




MP2/6-31+G*" 


2.616 


1.832 


1.083 


108.5 




B3LYP/6-311+G(2d,d,p) 


2.410 


1.931 


1.078 


105.5 




CCSD(T)/spd/g'' 


2.502 


1.853 


1.080 


107.6 


F/Br 


BH&HLYP/cc-pVTZ(+X) 


2.392 


2.052 


1.068 


105.2 




mPWH&HPW91/cc-pVTZ(+X) 


2.404 


2.007 


1.070 


106.4 




mPWlK/6-31+G* 


2.411 


2.022 


1.075 


106.2 




mP WIK /cc-pVDZ (+X) 


2.360 


2.045 


1.081 


105.2 




mPWlK/cc-pVTZ(+X) 


2.389 


2.022 


1.072 


106.0 




MP2/6-31+G*" 


2.528 


2.028 


1.081 


106.8 


CI/Br 


B3LYP/cc-pVTZ(+X) 


3.112 


2.024 


1.079 


106.7 




BH&HLYP/cc-pVTZ(+X) 


3.149 


1.986 


1.072 


107.3 




mPWlPW91 /cc-pVTZ(+X) 


3.088 


1.989 


1.080 


107.4 




mPWH&HPW91/cc-pVTZ(+X) 


3.113 


1.958 


1.074 


107.8 




B97/cc-pVDZ(+X) 


3.134 


2.021 


1.093 


106.8 




B97-l/cc-pVDZ(+X) 


3.098 


2.022 


1.093 


106.8 




HCTH/cc-pVDZ(+X) 


3.269 


2.012 


1.094 


107.1 




HCTH-120/cc-pVDZ(+X) 


3.132 


2.031 


1.093 


106.8 




mPWlK/6-31+G* 


3.070 


1.966 


1.079 


107.9 




mPWlK/cc-pVDZ(+X) 


3.114 


1.974 


1.085 


107.5 




mPWlK/cc-pVTZ(+X) 


3.109 


1.966 


1.076 


107.8 




MP2/6-31G* 


3.092 


1.992 


1.082 


107.6 




MP2/6-31+G*" 


3.199 


1.992 


1.084 


107.7 




CCSD(T) /sp<i/= 


3.095 


1.986 


1.082 


107.5 



« Prom Glukhovtsev, M.N.; Pross, A.; Radom, L.; J. Am. Chem. Soc. 1996, 118, 6273. 
^ From Schmatz, S.; Botschwina, P.; Stoll, H.; Int. J. Mass Spectrom. 2000, 201, 277. 

Prom Botschwina, P.; Horn, M.; Seeger, S.; Oswald, R.; Ber. Bunsenges. Phys. Chem. 
1997, 101, 387. 



lArJLiiii IV. (oreometnes ^^A , degree j oi 


ion-molecule complexes, 




j-V,Q G1\T 

oi tne k5iN2 


non- identity reactions (complex 3). 










Y/X Method/Basis Set 


r(X. • -C) 


r(C-Y) 


rfC-Hl 


ZHCY 




3.286 


1.425 


1.085 


108.8 


jDxi.06xi.i_f I x^ / lyC V X zj I 1 


3.268 


1.401 


1.077 


109.0 


mPWIPWQI /rr-nVTZC-l-X") 


3.246 


1.409 


1.085 


109.2 


mPWHJ^rHPWQl /rr-nVTZC+X") 

XXIX V ¥ XXOOXXX V V J X / Vw-V^, 1/ V X. ZJ l ^^./V 1 


3.230 


1.389 


1.079 


109.4 


RQ7 /rr nVD/Y-J-X") 


3.320 


1.432 


1.098 


108.7 


Rq7-1 /rr-nVnZl'+Xl 


3.268 


1.432 


1.098 


108.7 


HCTH/rr-nVDZl'+X") 

xxv.^ x_ XX / |j V iyzj I ^^yv y 


3.528 


1.431 


1.099 


108.7 


riw -L n- 1 cc-p V ljzj i^-|-j'v ^ 


3.358 


1.436 


1 nQQ 


1 flR 7 


mPWITC/R "^lO-n* 

xiir^ vv J--"^/ o-OJ-T^*^ 


3.241 


1.403 


x.uoo 


xuo.o 


TTiPWI K /rr-nVDZr-l-X^l 

XXXX V V XXV / l^L/ p V X-/ZJ I ^^./V 1 


3.252 


1.407 


1.089 


109.0 


mPWI K /rr- nVTZf^+X'l 

lUX V V XXV / k^l.;^ V X ZJl T^.^v 1 


3.235 


1.394 


1.081 


109.3 


iVLjT Z / 


3.227 


1.415 




1 HQ A 


MP2/6-31+G*" 


3.255 


1.438 


1.086 


108.0 


B3LYP /6-311G(2d,d,p) 


3.299 


1.418 


1.089 


109.4 


B3LYP /6-311+G{2d,d,p) 


3.271 


1.430 


1.087 


108.7 


CCST>n) / svdf 


3.188 


1.418 


1.086 


108.9 


F/Rr RST,YP/rr-nVT7C-4-X") 

X / ±JL -L_>OX_i XX/ (jl^ £J V J- ZJ I ^^./V 1 


3.497 


1.421 


1.085 


108.8 


RH;?rHT YP /rr-nVT7C-)-X"l 

XjXX06XXJ_( X X^ / V X Zj I J 


3.482 


1.398 


1.078 


109.0 


mPWIPWQI /rr-nVTZC-l-X") 

XXXX VVXX VVi7X / (jl^ p V J- ZJ I |^-/V 1 


3.451 


1.406 


1.086 


109.2 


mPWHX^HPWQI /rr-nVTZC+X") 

XXXX V V XXOoXXX V\ iJ ± / i^Kj V X. ZJ I ^^./V J 


3.433 


1.386 


1.080 


109.3 


RQ7 /rr nVD/Z+X") 


3.518 


1.428 


1.098 


108.7 


Rq7-1 /rr-nVnZl'-l-X'l 

xj t A- 1 \^\^ yj V x_/ ZJ ^ J V J 


3.458 


1.428 


1.098 


108.7 


HCTH /rr-nVnZf+X~) 

xxv.^ XXX/ ^^.j VJ y 1 ZJ I ^^./v / 


3.746 


1.428 


1.099 


108.8 


HGTH-1 90 /rr-nVnZI'-l-X'* 

XX^ X XX- XZU / p V X-ZZj ^ 


3.565 


1.433 


1.100 


108.7 


mx^ vv xxv/o-ox^^ 


3.367 


1.401 


X.UoO 


1UO.C5 


mP WIK /cc-pVDZ (+X) 


3.445 


1.404 


1.090 


109.0 


mPWlK/cc-pVTZ(-|-X) 


3.443 


1.391 


1.081 


109.3 


MP2/6-31G* 


3.315 


1.414 


1.088 


109.5 


MP2/6-31+G*'' 


3.457 


1.435 


1.087 


108.0 


Gl/Rr R3TYP/rr-nVTZ('-l-X'l 

Vw'i / X-JL xj oxj XX/ L^v_j Lj V J- ZJ ^ ^^yv J 


3.405 


1.839 


1.081 


108.0 


RHX^HT YP /rr-nVTZl'+X"; 

xjxxooxxxj 1. l / p V J- ZJ 1 »- / 


3.416 


1.813 


1.074 


108.3 


mPWIPWQI /rr-rtVT7,(4-X) 

XXXX VVXX vvi7x / (ji^ £y V J. ZJ I ^^yv J 


3.368 


1.812 


1.082 


108.5 


mPWHJ^^HPWQl /rr- nVTZC-l-X") 

XXXX V v XX OCX XX V v J X / l/ v x. zj i v i 


3.369 


1.791 


1.076 


108.7 


RQ7 /rr-nVr>7,(-i-X\ 
xjy ( / utj p V x-/zji T^./v I 


3.396 


1.850 


1.094 


107.8 


RQ7-1 /rr-nVnZC-l-X") 

xj£7 ( J. / i^v/ V x^ ZJ I ^^./v y 


3.354 


1.851 


1.094 


107.7 


HCTH/cc-pVDZ(+X) 


3.588 


1.832 


1.095 


108.2 


HCTH-120/cc-pVDZ(-|-X) 


3.440 


1.842 


1.095 


108.0 


mPWlK/6-31+G* 


3.274 


1.809 


1.080 


108.6 


mPWlK/cc-pVDZ(+X) 


3.381 


1.800 


1.087 


108.6 


mPWlK/cc-pVTZ(-|-X) 


3.371 


1.797 


1.077 


108.7 


MP2/6-31G* 


3.257 


1.812 


1.084 


109.1 



MP2/6-31+G*'' 3.457 1.807 1.085 108.9 

CrSD(T)/«/w//'' 3.318 1.820 1.083 108.2 

" From Glukhovtsev, M.N.; Pross, A.; Radom, L.; J. Am. Chem. Soc. 1996, 118, 6273. 
^ Prom Schmatz, S.; Botschwina, P.; StoU, H.; Int. J. Mass Spectrom. 2000, 201, 277. 

From Botschwina, P.; Horn, M.; Seeger, S.; Oswald, R.; Ber. Bunsenges. Phys. Chem. 
1997, 101, 387. 



lArJLiiii V. ^jeometries ^^A , aegreej oi tiie auiIsA 


transition structures, 


(^A=r , ana rsr ) 


of the SN2 identity reactions (complex 2). 






Species Method /Bosis Set 




rfG-Hl 

xxy 




1.808 


1.071 




1.854 


1.070 




1.823 


1.062 


mPWlPWQl /rr-nVTZ(^+X1 


1.824 


1.071 


mPWHX^HPWQI /rr-nVTZC+X'l 


1.797 


1.064 




1.852 


1.085 


B97-l/cc-pVDZ(+X) 


1.848 


1.084 


HriTH / rr-nVr)7.( 4-X) 

XX v-' X. XX/ ^ji^^ V zj 1 ^^j^v ; 


1.875 


1.085 




1.872 


1.086 


mPWI K /fi-'il 4-n* 


1.807 


1.070 




1.810 


1.077 


mx^ w ixY/ cc-p V i Zi ^_"r-A. 


i.c5U4 


l.UoD 


MP2/6-31G* 


1.778 


1.076 


MP2/6-31+G*" 


1.837 


1.074 


B3LYP/6-311G(2d,d,p) 


1.830 


1.074 


B3LYP /6-311+G(2d,d,p) 


1.862 


1.073 


X— Gl GGSn(^TWrr-r)VOZ4-1 


2.305 


1.070 


R3T YP /rr-nVTZ(^4-X'l 


2.355 


1.069 


jjn.o6n.xj X r / |J V J.zjiT^y\.i 


2.332 


1.061 


mPWlPWQl /rr-nVTZr+X'l 

iixx vvxx vvcx/*^^ V X Z—J I 1 J V J 


2.310 


1.070 


mPWH^HPWqi /rr-nVTZC+X'l 

lllX VVXX Oo XXX y \ Ij A. 1 K^V^ \J V X ZJ I ^^./V 1 


2.290 


1.064 


RQ7 /rr- nVr)7('4-X'l 


2.349 


1.083 


RQT-l /cT-nVDZf+X~l 

xj ci x/i.^^^ sr \ 1 


2.344 


1.083 


HGTH /rr-nVDZ (+X\ 

XX v.^ X XX / VjK^ \j V x^ ZJ I ^^^v 1 


2.374 


1.084 


XX\_i X XX- XZU / p V X_/Zj1 T^TV I 


2.372 


1.084 


lllX^ VVXXV/U OX^VJ 


2.313 


1.069 


lux V V xxv / i^k^ jj V x_/iji T^yv 1 


2.303 


1.076 


Illr VV ilv/ CC-p V IZjl-f-A.) 




1 flAC^ 
l.UDO 


MP2/6-31G* 


2.308 


1.072 


MP2/6-31+G*'^ 


2.317 


1.073 


B3LYP/6-311G(2d,d,p) 


2.366 


1.072 


R3T YP /fi-31 1 -ir-Ci(1(\ H n'l 

xjtjxj 1 L / (J ijx x^^vjr 1 ^u^Uj^J 1 


2.364 


1.072 


X— Rr GGSn(^TWrr-r)VOZ4-1 

-TV XJX Ky\JvJ±-/ I X. 1 / I^L/ ^ V Vi^ZJ^^ X 


2.461 


1.071 


R3LYP /rr-nVTZr+X'l 

J-J<JXJ XX/ W V X ZJl ^^^V J 


2.511 


1.069 


RHJ^^HLYP /rr-nVTZr+X'l 

J— 'XXOOXXX-J XX/ U V X ZJl ^^^V J 


2.488 


1.062 


mPWlPW91 /cc-pVTZ(+X) 


2.464 


1.070 


mP WH&HPW91/cc-pVTZ(+Xj 


2.442 


1.064 


B97/cc-pVDZ(+X) 


2.502 


1.084 


B97-l/cc-pVDZ(+X) 


2.496 


1.084 


HGTH/cc-pVDZ(-|-X) 


2.533 


1.085 


HGTH-120/cc-pVDZ(+X) 


2.530 


1.085 



mPWlK/6-31+G* 2.430 1.069 

mPWlK/cc-pVDZ(+X) 2.459 1.076 

mPWlK/cc-pVTZ(H-X) 2.447 1.066 

MP2/6-31G* 2.444 1.073 

MP2/6-31+G*" 2.480 1.074 

« From Glukhovtsev, M.N.; Pross, A.; Radom, L.; J. Am. Chem. Soc. 1995, 117, 2024. 



lArJLiiii Vi. ^jeometries ^^A , aeg 


^ree ) oi tne i Uri3A 


transition 


structures, 


+ViQ Q1\T 

OI tne k5i\2 


non- identity reactions (complex 2). 










Y/X Method/Basis Set 


t(Y- ■ -C) 






ZHCX 




2.143 


2.091 


1.070 


99.0 




2.065 


2.107 


1.062 


97.2 


mPWI PWQl /rr-nVTZf^+X'l 


2.072 


2.086 


1.071 


97.8 


mPWHJl^HPWqi /rr-nVTZr+X'l 


2.016 


2.088 


1.064 


96.7 


RQ7 /rr nVn7('-l-X'> 


2.143 


2.091 


1.085 


98.8 


RQ7-1 /rr-nVDZr+X") 


2.139 


2.086 


1.085 


98.8 


HCTH /rr-nVnZl'+X'l 

J-AX^* 1. 1.1. 1 KjK^^ V 1.J Li \ J 


2.190 


2.098 


1.086 


99.2 


HflTH-1 9n/rr-nVnZ('+X~l 


2.208 


2.077 


1.086 


100.0 


mPWl K /fi-'il 4-0* 

lUA VV l.L\./ yj O -L 1^ vjf 


2.063 


2.080 


1.070 


98.0 


nil VV -L-iV/ CC-p V LJZjy-\-J^ ) 




Z.UC50 


L.U to 


Q7 A 


mPWlK /cc-pVTZ (+X) 


2.029 


2.089 


1.066 


97.0 


MP2/fi-S1-i-0*" 


2.016 


2.142 


1.073 


95.6 


RST/YP /6-Sl 1 +Gf2d d n) 

i—f 'J J_J ± ± / ^ J- -L V_T I ^^4. 1 


2.173 


2.080 


1.074 


99.7 


CCSDl'Tl/sodfo'' 


2.030 


2.121 


1.072 


96.3 


P/Rr RHXrHTYP /rr nVTZl'+X") 


2.174 


2.181 


1.063 


100.0 


mx^ VV xio6xir^ VV y 1/ CC-p V LLiy-T^) 


z.iuo 


9 1 7c: 


1 C\fiK 

l.UOO 


QQ 9 

yy.z 


mPWlK/6-31+G* 


2.114 


2.181 


1.070 


99.2 


mPWlK/cc-pVDZ(+X) 


2.179 


2.145 


1.078 


100.9 


mPWlK /cc-pVTZ (+X) 


2.129 


2.168 


1.067 


99.7 


MP2/6-31+G*" 


2.108 


2.242 


1.075 


97.9 


ni/Rr RST,YP/rr-r)VTZf-t-X~) 


2.416 


2.451 


1.069 


91.8 


RRX^RTYP /rr nVTZl'+X") 


2.388 


2.433 


1.062 


91.6 


mPWIPWQI /rr-nVTZf+X") 


2.369 


2.409 


1.070 


91.7 


mPWHJl'HPWQl /rr-nVTZf+X"! 

xxxx V V xxcoxxx V V t/ x / v i zj \ v j 


2.342 


2.391 


1.064 


91.5 


RQT/rr-nVnZC-l-X") 


2.408 


2.445 


1.084 


91.7 


RQ7-1 /rr-nVnZC+X'l 


2.403 


2.439 


1.084 


91.7 


HCTH /rr-nVDZr+X") 


2.488 


2.467 


1.085 


91.3 


WPTTT 19n/<-*- T->Vr>7/'-U"X"'l 
n v_/ _L n- ±zu / cc-p v uzji t^yv i 






J-.UOO 


yi.o 


lilJr VV llv / U-O J-^*^ 


Z.OIO 


2 421 


-L .uoy 


90 1 


mPWlK/cc-pVDZ (+X) 


2.363 


2.401 


1.076 


91.8 


mPWlK/cc-pVTZ(+X) 


2.349 


2.396 


1.065 


91.6 


MP2/6-31G* 


2.336 


2.419 


1.072 


90.5 


MP2/6-31+G*" 


2.371 


2.430 


1.073 


91.4 


CCSD(T)/spd/= 


2.354 


2.422 


1.072 


91.2 



« Prom Glukhovtsev, M.N.; Pross, A.; Radom, L.; J. Am. Chem. Soc. 1996, 118, 6273. 
^ From Schmatz, S.; Botschwina, P.; Stoll, H.; Int. J. Mass Spectrom. 2000, 201, 277. 

Prom Botschwina, P.; Horn, M.; Seeger, S.; Oswald, R.; Ber. Bunsenges. Phys. Chem. 
1997, 101, 387. 



